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of the qubits in a quantum computer, one can find eigenvectors 
very quickly. If a unitary operation U is implemented as a time- 
independent Hamiltonian, for instance, one can collapse the state 
of the computer to a nearby eigenvector of U with a measurement 
^ ! of the energy. We examine some recent proposals for quantum 

\ computation using time-independent Hamiltonians and show how 

to convert them into "artificial orbitals" whose energy eigenstates 
q ' match those of U. This system can be used to find eigenvectors 

and eigenvalues with a single measurement. We apply this tech- 
q . nique to Grover's algorithm and the continuous variant proposed 

by Farhi and Gutmann. 
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1 Abrams and Lloyd's Eigenvector Algorithm 

When Feynman first suggested quantum computers, he envisioned them as 
being used to simulate actual physical systems efficiently Since the size 
of the Hilbert space increases exponentially with the number of particles, 
problems quickly become intractible for classical computers. Lloyd showed 
[9] that most physical systems can be simulated in polynomial time with the 
product of few-qubit unitary matrices. Abrams and Lloyd later presented a 
quantum algorithm [1] for finding eigenvectors and appplied it to physical 
systems. The rate of convergence, however, depends on the eigenvalue; for 
eigenvectors whose eigenvalues are exponentially close to a multiple of 2tt, 
even a relatively good guess will take an exponentially long time to converge. 

Their algorithm attempts to find eigenvectors of some unitary matrix U 

k 

by taking a guessed eigenvector > and calculating J2 U^\ip >. If \ip > is 
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1 



near an eigenvector \<f> > of U, that eigenvector will dominate in the sum. 
The rate of convergence depends on three factors: the initial guess, k, and 
the eigenvalue of \<f>>. Consider, for example, a very small rotation about 
an axis. Each application of U only moves the vector by a small amount; the 
sum will stay very close to the original guess for small k. In order to find 
the axis exactly, one would need k large enough that ip gets rotated all the 
way around the axis. For exponentially small eigenvalues, that would make 
k exponentially large. 

2 Artificial Orbitals 

There is another way to find eigenvectors. The eigenvectors of any unitary 
operation U = Um ■ ■ ■ U 2 Ui are also those of some measurement. If one re- 
stricts measurements of a quantum computer to measurements of the logical 
qubits, there is no simple relation between the eigenvectors of U 2 and U\ and 
their product, so transforming between an initial basis and the U basis is not 
easy. Allowing measurements of other observables, however, can increase the 
power of the computer. If U is implemented as a time-independent Hamilto- 
nian, for example, a measurement of the energy of the system will collapse 
the state onto a nearby eigenvector of U . 

Standard models of quantum computation apply different Hamiltonians 
for each step of the computation. The Hamiltonian describing the entire evo- 
lution of the system is therefore highly time dependent. Recently, however, 
two different groups have proposed a new paradigm for quantum compu- 
tation [7,8,10,11]. Their model is very much like a complete quantum net 
[12]; instead of applying the operations U t to a quantum state in a single 
position over time, as in NMR quantum computation [3] and ion traps [2], 
the quantum state moves through the net, which is static. Annihilation and 
creation operators move the state from one part of the net to another and 
are interspersed with few-qubit interactions U x . 

The system can be viewed as having two main registers: the "virtual time 
step" register, which tells how far through the quantum net the computation 
has progressed; and the qubit register, in which the computation occurs. If P x 
projects the system onto virtual time step x, one can view the Hamiltonian 
as implementing the system of equations 

p x +i\il»= u x+1 p x \ij> 
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for each step x in the quantum net. 



The physical implementations suggested by these two groups replace the 
virtual time step register with the position of the electrons in a quantum 
circuit. We can adapt the model, however, to physical implementations 
where the quantum state does not move through space. For implement- 
ing an M-step algorithm, one can imagine a system composed of the sum 
of M conditional Hamiltonians plus a "stepping" Hamiltonian that are all 
time-independent. The conditional Hamiltonians H t only change the quan- 
tum state when the virtual time step register is in the state t. They act 
on the lg{M) + k qubits comprising the virtual time step register plus k for 
implementing Ut, the operation at that time step; since virtually any two- 
qubit gate is universal, we can take k = 2. The stepping Hamiltonian simply 
increments the time step, acting on lg{M) qubits. Once the virtual time step 
has reached M, one measures the qubit register and reads the output. 

If we allow the stepping Hamiltonian to loop back to zero after reaching 
M, we make the output of the circuit feed back into the input. The quantum 
net or circuit becomes a cylinder, and the energy eigenstates of the system 
take the form 



where is the jth cyclic permutation of U — Um ■ ■ ■ U 2 U±, and ipui is an 
eigenvector of UK 

We can use such a system to find eigenvalues and eigenvectors of U. First, 
we initialize the state with our best guess of an eigenvector. Measuring the 
energy of this system will collapse the state onto a nearby eigenvector and 
give us the corresponding eigenvalue; measuring the virtual time step register 
next collapses the state to a random time step. We repeat these three steps 
until we get virtual time step zero; since there are M time steps, this will take 
O(M) repetitions. At this point, the qubit register contains an eigenvector 
of U. The quality of our initial guess and M are the only logical parameters 
governing the success rate; it is independent of the eigenvalue. 
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3 Eigenvectors of Quantum Search Algorithms 



Because of its eigenvalue independence, one particularly exciting applica- 
tion of this process is finding eigenvectors of quantum search algorithms like 
Grover's algorithm [6] and Farhi and Gutmann's algorithm for the continuous 
Grover problem [4]. 

Grover's algorithm is usually written G = I S I W = Ih\q>Iw — 

HI HI W . In 

this notation, I k is an inversion that negates the \k> component of the wave 
function, and H is the Walsh-Hadamard transform. I w is usually referred to 
as the "quantum oracle," since it "marks" the unknown state \w >, while I s is 
known as the "amplitude amplification," because it increases the probability 
of measuring a marked state. If the quantum computer begins in the state 
|s>, 0{s/N) applications of G will rotate the state into \w>. 

Farhi and Gutmann replace the quantum oracle with a time-independent 
Hamiltonian having a single excited state \w > and the amplitude amplifier 
with a similar Hamiltonian having a single excited state \s>. In this system, 
the quantum state will evolve from \s> to |iu> in 0(y/~N) time. 

These algorithms deal with three different bases. One is the logical qubit 
basis in which the oracle's test for the marked state occurs. If the oracle 
implmented a test for, say, the satisfiability of a given circuit, the logical 
qubit basis would have each qubit represent one input to the circuit, and 
w is that assignment of inputs that satisfies the circuit. The state \w > is 
an eigenvector of the logical qubit basis. The second basis is the Hadamard 
transform of the logical qubit basis. This basis is conjugate to the first in the 
same way that momentum is conjugate to position, or spin measured along 
the x axis is conjugate to spin measured along the z axis. The state \s> is 
an eigenvector of this conjugate basis. 

The third basis is the "search" basis; that is, the vectors comprising this 
basis are those vectors that only change by a global phase when the search 
algorithm is applied. Geometrically, Grover's algorithm G is a rotation in the 
(\s>, \ w>) plane, while Farhi and Gutmann's algorithm F is the generator 
of a reflection in the same plane. Eigenvectors of G that span the plane are 
the circularly polarized states 

1 „ , 
\±G >— —i=(\s> ±i\w>), 



4 



while eigenvectors of F that span the plane are diagonally polarized: 

\± F >~ ^(|s> ±k>). 

The initial state \s> is a superposition of |±> that constructively inter- 
feres in the \s> direction: 

1 n 

l«>- ^fd +> +!->)• 

Applying the search algorithm rotates the phase of these two eigenvectors in 
opposite directions, gradually shifting the direction of constructive interfer- 
ence to the \w> direction: 

1 n 

F>- -7=0 + > -|->)- 



All the analyses of the optimality of Grover's algorithm (see, for example, 
the references in [5]) show that the phase rotation cannot happen faster than 
0(y/~N). This is true. However, since \s> is a superposition of two eigenvec- 
tors of the search basis, a measurement in the search basis will collapse \s> 
onto one of them. Instead of waiting an exponentially long time for the state 
to evolve from \s> to \w>, we should measure in the search basis instead. 
This will instantly collapse the wave function onto one of |± > with certainty. 

Since F is assumed to be implemented with time-independent Hamiltoni- 
ans, its energy basis matches the search basis. A measurement of the energy 
will collapse the state \s> onto one of |ztp> with certainty. A subsequent 
measurement in the logical qubit basis will yeild \w> with probability 1/2. 
Had we left the system alone, the system would have eventually arrived at 
the state \w>, it is true, but this way we can get there with only two mea- 
surements half of the time. 

Farhi and Gutmann did not suggest an implementation for F; however, 
we can implement G with a small quantum circuit as an artificial orbital 
and solve the problem in much the same way. One loop around the artifi- 
cial orbital will change the quantum state by one iteration of G. Since the 
Hamiltonian is time-independent, there is an implicit solution to the problem 



5 



in the energy basis. We set up the initial state \s> and measure the energy; 
this collapses the internal state of the computer to a superposition \xjj> that 
satisfies Pi\ip>= \t%> |±g>- Next, we measure the virtual time step. If we 
do not get t%, then we repeat the previous three steps. The last measurement 
succeeds with probability 1/M. Finally, we measure the qubits; since |±<j> 
is equal parts \s> and \w >, with probability 1/2 we find the qubits in the 
state \w>. Thus we have a search whose success probability depends on the 
number of bits only insofar as it affects the size of the circuit. For Grover's 
algorithm, the size is at most a polynomial in the number of bits. 



Conclusion 

By allowing measurements of observables other than the values of the qubits, 
one can find eigenvalues of arbitrary products of M few-qubit unitary matri- 
ces in O(M) time. This includes simulations of most physical systems as well 
as quantum search algorithms, allowing one to solve NP-complete problems 
in polynomial time. 
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